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The Navier-Stokes transport coefficients of a granuiar dense fluid of smooth ineiastic hard disks 
' or spheres are expiicitiy determined by soiving the inelastic Enskog equation by means of Grad's 

moment method. The transport coefficients are explicitly determined as functions of the (constant) 
■ coefficient of restitution and the solid volume fraction. In addition, the cooling rate is also calculated 

£\j ' to first order in the spatial gradients. The calculations are performed for an arbitrary number of 

dimensions. The results are not limited to small dissipation and are expected to apply at moderate 
densities. It is found that the expressions of the Navier-Stokes transport coefficients and the cooling 
rate agree with those previously obtained from the Chapman-Enskog method by using the leading 
terms in a Sonine polynomial expansion. This shows the equivalence between both methods for 
f"^ . granular fluids in the Navier-Stokes approximation. A comparison with previous results derived 

£\j ■ from Grad's moment method for inelastic disks and spheres is also carried out. 
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I. INTRODUCTION 



It is well established that granular media under rapid flow conditions admit a hydrodynamic description. 1 The 
essential difference from that for ordinary fluids is the absence of energy conservation, yielding subtle modifications 
of the conventional Navier-Stokes (NS) equations for states with small spatial gradients of the hydrodynamic fields. 
To gain some insight into the general description of the dynamics of grains, a simple model is usually considered for a 
granular gas: a system constituted by smooth hard spheres or disks with inelastic collisions. The loss of energy in each 
binary collision is accounted for by a constant coefficient of normal restitution a < 1, the case a — 1 corresponding 
to elastic collisions. 

At a kinetic theory level, the Boltzmann 2 and Enskog 3, equations conveniently adapted to account for inelastic 
binary collisions, have been employed in the past few years as the starting points to derive the NS hydrodynamic 
equations from a more fundamental point of view. In particular, assuming the existence of a normal or hydrodynamic 
solution 5 for sufficiently long space and time scales, the Chapman-Enskog method 6 has been applied to calculate 
the distribution function /(r, v,t) through first order in the spatial gradients. Use of this distribution allows one to 
determine the dependence of the NS transport coefficients on the coefficient of restitution for dilute 7 and moderately 
dense 8,9 gases. In contrast to previous attempts, 10-12 the results derived in Refs. 7-9 do not impose any constraint on 
the degree of dissipation and take into account the (complete) nonlinear dependence of the transport coefficients on a. 
However, as for elastic collisions, 6 the exact forms of the NS transport coefficients require the solution of a set of linear 
integral equations and so the leading terms in a Sonine polynomial expansion (first Sonine approximation) are usually 
considered to get explicit expressions for all the above coefficients. In spite of this simple approach, the corresponding 
' analytical results compare quite well with Monte Carlo simulations, 13 except at high dissipation for the heat flux 
, transport coefficients. 14,15 Motivated by this disagreement a modified version of the first Sonine approximation has 
£N) ■ been recently proposed 16,17 that significantly improves the a-dependence of the heat flux transport coefficients and 
t-H | corrects the discrepancies between simulation and theory. 

An alternative procedure to solve the Boltzmann equation for a dilute gas is by means of the moment method. 
The objective of the method is to evaluate the velocity moments of the distribution function /(r, v, t) rather than its 
explicit form as in the Chapman-Enskog method. Those velocity moments provide an indirect information on / and, 
additionally, its first few moments (the mass density p, the flow velocity U, the temperature T, the pressure tensor 
Pij, and the heat flux vector q) characterize the macroscopic state of the gas. Nevertheless, when one multiplies 
both sides of the Boltzman equation by a given set of functions ^"(v) and integrates over the molecular velocity, in 
general one obtains an infinite hierarchy of moment equations. 5 This infinite hierarchy can be recursively solved for 
some specific interactions potentials (such as Maxwell molecules, namely, when particles repel each other with a force 
inversely proportional to the fifth-power of the distance) and/or some special non-equilibrium situations. 18 On the 
other hand, beyond this simple interaction potential, one has to resort to an approximate method to solve the above 
hierarchy of moment equations. The most widely used method was devised by Grad more than fifty years ago. 19 
The idea of Grad's method is to assume /(r, v, t) to be a local Maxwellian /m(i", v, t) times a sum over Hermite 
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polynomials i?fc(v), i.e., 

N-l 

f(r,v,t)^ f M (r,v,t)J2 Cfc(r,t)H fc (v). (1) 

fc=0 

The coefficients appearing in each of the velocity polynomials (v) are chosen by requiring that the corresponding 
velocity moments of Grad's solution (1) be the same as those of the exact velocity distribution function. There are 
N arbitrary quantities which may be identified with the basic (p, U, T, Pij, and q) and higher moments and can be 
determined by recursively solving the corresponding transfer equations for the above moments. A reasonable choice 
for a three-dimensional ordinary dilute gas is N — 13; in such a case the unknowns are the basic hydrodynamic fields 
(p, U, and T) and the irreversible momentum and heat fluxes (Pij —nTSij and q). In this case, the method is referred 
to as Grad's thirteen moment method. 5 In the case of a general dimensionality d the number of involved moments 
is d(d + 5)/2 + 1. Although Grad's moment method was originally proposed for dilute gases, its extension to dense 
gases is easy since one only has to consider the kinetic contributions to the fluxes in the trial solution (1). 

Grad's moment method has been also applied to granular gases. In the context of the inelastic Enskog equation, 
Grad's thirteen-moment method was employed several years ago by Jenkins and Richman 20 ' 21 to determine the stress 
tensor, the heat flux, and the cooling rate in the NS approximation (linear theory). Although the application of 
Grad's method to the Enskog kinetic equation is not restricted to nearly clastic particles (a ~ 1), the results derived 
by Jenkins and Richman 20 ' 21 neglect the cooling effects on the granular temperature T due to the cooling rate. Given 
that this assumption can only be justified for nearly elastic systems, their expressions for the NS transport coefficients 
differ from those obtained in Refs. 8,9 from the Chapman-Enskog method for arbitrary degree of inelasticity. More 
recently, 22 Grad's method has been also applied for weakly inelastic dilute gases with a coefficient of restitution which 
depends on the relative velocity. 23 

The aim of this paper is to use Grad's method to obtain the NS transport coefficients of <i-dimensional granular 
fluids described by the Enskog kinetic equation. This study extends a previous work of the author 24 for dilute 
granular gases and so it provides a description of hydrodynamics and transport at higher densities. With respect to 
the Jenkins-Richman theory, 20 ' 21 the present results incorporate two new ingredients not considered in the previous 
study. First, they are obtained by accounting for the time dependence of the temperature coming from the inelastic 
cooling. As a consequence, the corresponding expressions of the NS transport coefficients hold for arbitrary degree of 
dissipation. Second, a new scalar field (the full contracted moment of fourth order c) is added to the usual thirteen 
moments of mass density, velocity, temperature, and the kinetic contributions to the pressure tensor and heat flux 
vector. The inclusion of the moment c in Grad's solution has been already considered in previous works 25,26 on dilute 
granular gases. A direct consequence of the presence of the field c is a new contribution to the cooling rate proportional 
to the divergence of the flow velocity as well as several new contributions to the transport coefficients coming from 
non-Gaussian corrections to the distribution function in the homogeneous cooling problem. 27 The question arises then 
as to wether, and if so to what extent, the conclusions drawn by Jenkins and Richman 20 ' 21 may be altered when the 
above two new ingredients are incorporated in Grad's solution. 

In the NS approximation (namely, when only linear terms in the spatial gradients are retained), the constitutive 



equations for the pressure tensor P^, the heat flux q and the cooling rate C are 3 ' 8 ' 20,21 

Pij = pS tJ - v [diUj + djU, - ^SijV • u) - 7 <%V • U, (2) 

q = -kVT - pVn, (3) 

C = Co + CuV • U, (4) 



where p is the hydrostatic pressure and n is the number density. In addition, r\ is the shear viscosity, 7 is the bulk 
viscosity, k is the thermal conductivity, p is a new transport coefficient not present for elastic collisions, and Co and 
Cu are the zcroth- and first-order contributions, respectively, to the cooling rate. The above NS transport coefficients 
can be written in reduced forms as 

p{a,4>) =nTp*(a,(t>), r)(a, 0) = r?(l, 4>)rf{a, <f>), j(a, <j>) = r?(l, 0)7* (a, 4>), (5) 



K(a,4>) = k(1,4>)k* (a, <f>), p(a,4>) = ■ — p*(a,4>), Co (<*,<£) = n tt Co (a,4>), (6) 

n T7(l, 0) 
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whore 

TT d / 2 

-na d , (7) 



2 d - 1 dr(f)' 

is the solid volume fraction and a is the diameter of the spheres. The coefficients 7/(1, (f>) and <j>) are the values of 
the elastic shear viscosity and thermal conductivity, respectively, given by the Enskog equation. 28 The results derived 
in this paper show that the dimcnsionless coefficients p* , n* , 7*, k*, /U*, Q, and (u are in general nonlinear functions 
of the coefficient of restitution a and the solid volume fraction (p. In addition, the expressions of all the above 
dimensionless NS transport coefficients obtained here from Grad's moment method agree with those derived from 
the Chapman-Enskog expansion in the first Sonine approximation. 8,9 This confirms the expected mutual consistency 
between both methods for solving the inelastic Enskog equation in the NS domain. 

The plan of the paper is as follows. In Sec. II the Enskog kinetic equation and associated macroscopic conservation 
laws for a granular fluid are introduced. An overview of Grad's moment method used for solving this kinetic equation 
is given in Sec. III. The explicit results for the NS transport coefficients are provided in Sec. IV, with the details of the 
calculations appearing in Appendices A, B, and C. A comparison with the results obtained by Jenkins and Richman 
for disks 20 and spheres 21 is done in Sec. V, showing significant discrepancies between both theories especially in the 
heat flux transport coefficients. Finally, a short summary of the results derived in the paper is presented in Sec. VI. 

II. ENSKOG KINETIC THEORY AND CONSERVATION LAWS 

We consider a granular fluid composed by smooth inelastic disks or spheres of mass to and diameter a. The 
inelasticity of collisions among all pairs is accounted for by a constant coefficient of normal restitution < a < 1 
that only affects to the translation degrees of freedom of grains. The particular value a — 1 corresponds to elastic 
collisions (ordinary fluids). At a kinetic theory level, all the relevant information on the state of the system is given 
by the one-particle velocity distribution function /(r, v, t). For moderate densities, the inelastic Enskog theory 4 gives 
the time evolution of f(r,v,t). In the absence of an external force, it has the form 

(ft+vV)/(r,v,t) = J E [r,v|/(t)], (8) 

where 

J E {r,v 1 \f(t)]=o- d - 1 J dv 2 J d^e(a- g )0-g)[a- 2 f^(r,r- < 7,v[,V 2 ;t)-f^(r,r + ( T,v 1 ,v 2 -t)'\ (9) 

is the Enskog collision operator. In Eq. (9), d is the dimensionality of the system (d — 2 for disks and d = 3 for 
spheres), er = crer, er being a unit vector along the centers of the two colliding spheres, 8 is the Hcaviside step function, 
g = vi — v 2 is the relative velocity and 

(n , r 2 , vi , v 2 , t) = x(ri , r 2 |n(t))/(n , vi , t)/(r 2 , v 2 , t). (10) 

The primes on the velocities in Eq. (9) denote the initial values {v'i,v 2 } that lead to {vi,v 2 } following a binary 
collision: 

v' 1 =vi-i(l + a- 1 )(?-g)?, v 2 = v 2 + ^ (l + a" 1 ) (£ • g)c?. (11) 

Furthermore, x[r, r + cr|n(i;)] is the equilibrium pair correlation function at contact as a functional of the non- 
equilibrium density field n(r, t) defined by 

n(r,t) = J dv/(r,v,t). (12) 
The first d + 2 velocity moments of /(r, v, t) define the number density n(r, t), the flow velocity 

UM) = ^/ dvv/fcv,*), 

and the granular temperature 



(13) 
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where V(r, t) = v U(r, t) is the peculiar velocity. 

The exact macroscopic balance equations for n(r, t), U(r, t) and T(r,t) follow directly from the Enskog equation 
(8) by multiplying with 1, mv, and \mv 2 and integrating over v. After some algebra, one gets 8 

An + nV-U = 0, (15) 
pDtUi + djPij = 0, (16) 

AT + J- {dm + PijdjUi) = -QT, (17) 

where D t = dt + U • V is the material derivative and p = mn is the mass density. The cooling rate ( is (essentially) 
proportional to 1 — a 2 and is due to dissipative collisions. The pressure tensor P(r,i) and the heat flux q(r,i) have 
both kinetic and collisional transfer contributions, i.e., P = P fe + P c and q = q fc + q c . The kinetic contributions are 
given by 

P k (r,t) = J dvmW/(r,v,{), q fe (r,t) = | ^y 2 V/(r,v,i), (18) 



(19) 



and the collisional transfer contributions are 8 

P c (r, t) = ^-^-m<r d J dv 1 J dv 2 J da Q{a ■ g)(a ■ gfaa J dx f (2) [r - xa, r + (1 - x)a, vi, v 2 ; t] , 

q c (r,t) = ^-^ma d j dw 1 J dv 2 j da &{a ■ g)(a • g) 2 (G • a)a J dx / (2) [r - xa, r + (1 - x)a, vi, v 2 ; t] , (20) 

where /( 2 ) is defined in Eq. (10) and G = \ (V\ + V 2 ) is the velocity of the center of mass. Finally, the cooling rate 
is given by 

C(r,*)= ^^t 771 ^ j dVl J dv2 J d9e{a •g)( ? -g) 3 / (2) (r,r + < T,v 1 ,v 2 ;t). (21) 

Apart from the balance equations for the hydrodynamic fields n, U, and T, to completely characterize the macro- 
scopic state of the granular fluid one should also to derive the corresponding balance equations for the pressure 
tensor and the heat flux. These equations will be derived in the next section by considering an explicit form for the 
distribution function f(r,v,t). 

III. GRAD'S MOMENT METHOD. LINEAR THEORY 

Needless to say, to close the balance hydrodynamic equations (15)— (17) one needs to know the functional dependence 
of the momentum and heat fluxes and the cooling rate on the hydrodynamic fields n, U and T. A possible way of 
obtaining this dependence is solving the Enskog equation by means of the Chapman-Enskog method. 6 This was the 
procedure followed in Refs. 8,9 to determine Py, q and ( to first order in the spatial gradients (NS hydrodynamic 
order). Here, a different procedure will be followed: Grad's method of moments. 19 Although the method was originally 
devised to solve the Boltzmann equation for dilute molecular gases, here it will be used to determine the NS transport 
coefficients of a granular dense fluid described by the inelastic Enskog equation (8). 

As mentioned in the Introduction, Grad's moment method is based on the expansion of the velocity distribution 
function in a complete set of orthogonal polynomials (generalized Hermite polynomials), the coefficients being the 
corresponding velocity moments. However, given that the (infinite) hierarchy of moment equations is not a closed set 
of equations, one has to truncate the above expansion after a certain order. After this truncation, the above hierarchy 
of moment equations becomes a closed set of coupled equations which can be solved. This allows one, for instance, to 
get the explicit forms of the NS transport coefficients. An interesting question is to assess the differences between the 
results derived from the Chapman-Enskog expansion and Grad's moment method when only terms up to first order 
in the spatial gradients are retained in the constitutive equations for Py, q, and (. 

In the application of the standard Grad moment method for a dense fluid, the retained moments are the hydrody- 
namic fields (n, U, and T) plus the kinetic contributions to the irreversible momentum and heat fluxes (P^ — nTSij 
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and q fc ). In the three-dimensional case (d — 3), this implies that there are 13 moments involved in the form of the 
velocity distribution function /; hence this method is referred to as the 13-momcnt method. On the other hand, 
since we are interested in comparing the present results with those obtained for granular dense gases 8,9 from the 
Chapman-Enskog method, the full contracted moment of fourth order 



d(d + 2) 



AnT 2 



dv V\f(V) d(d + 2) 



(22) 



will be also included. The inclusion of the scalar field c to the thirteen moments of mass density, velocity, pressure 
tensor, and heat flux vector will allow us to make a close comparison with the previous Chapman-Enskog expressions 
derived for the cooling rate and the NS transport coefficients. 8 ' 9 

Under the above conditions, the explicit form of the non-equilibrium distribution function /(r, v, t) can be written 

as 



/(V) -+ / M (V) 

where 



/m(V,=„(^) J '%— V- (24) 
is the local equilibrium distribution function, 



(23) 



mV 2 d + 2\ fmV 2 \ 2 d + 2 mV 2 d(d + 2) 



s(v) = (^-— v, e{v) = 1 — ) — — — + ^-r^, (25) 



and 



IF, - P* - nTSij (26) 

is the traceless part of the kinetic contribution to the pressure tensor. The coefficients appearing in each one of the 
velocity polynomials in Eq. (23) have been chosen by requiring that the basic hydrodynamics fields (n, U, and T), 
the kinetic contributions to the pressure tensor and the heat flux vector, as well as the fourth contracted moment c 
of the trial function (23) to be the same as those for the exact velocity distribution function /. In addition, based 
on previous results derived in the NS order for dense gases, 8 it is expected that in the linear theory the distribution 
function (23) must depend on the divergence of flow velocity through the moments acting as coefficients. Given that 
the tensor IF, is traceless, it is expected that IF, is proportional to the tensor diUj +djUi — %5ij V • U while the heat 
flux q must be proportional to the vectors Vn and VT. Consequently, the only contribution proportional to V • U 
in Grad's solution (23) comes from the contracted fourth moment c. This necessarily implies (as we will show later) 
that the coefficient c is decomposed as 

c = c + c 1 u~ 1 W -V, (27) 

where v is an effective collision frequency for hard spheres [see Eq. (34) below]. According to the decomposition (27), 
while the coefficient c characterizes the deviations of / from its Gaussian form in the homogenous cooling state, 
the coefficient c\ gives rise to a new first-order contribution to the cooling rate £. This new contribution was not 
accounted for in previous results 20,21 based on Grad's method. 

For clastic collisions (a = 1), the coefficients Co and c\ vanish and so, one recovers the conventional form for the 
trial distribution / in Grad's thirteen moment method. 19 Note that, for the sake of simplicity, in the fourteen-moment 
approximation (23) some third-degree moments not included in the kinetic heat flux q fc have been left out. 25,26 The 
same can be said of the remaining polynomials of fourth order. The inclusion of the above moments would modify, 
for instance, the form of the cooling rate. However, as mentioned before, since we are interested in a theory with the 
same degree of accuracy as the one reported 8,9 by using the Chapman-Enskog method, only the kinetic contributions 
to the pressure tensor and the heat flux, as well as the contracted fourth order moment (22), will be included in the 
Grad's distribution function (23). 



IV. NAVIER-STOKES TRANSPORT COEFFICIENTS 



The form of the constitutive equations for the momentum and heat fluxes and the cooling rate in the NS order are 
given by Eqs. (2)-(4), respectively. In this section, the NS transport coefficients and the cooling rate will be explicitly 
determined by using Grad's distribution (23). For the sake of clarity, most of the technical details involved in these 
calculations are relegated to Appendixes A, B, and C and only the final expressions for rj, 7, k, /1, Co, and Cu will be 
displayed here. Let us consider first the momentum and heat fluxes. 
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A. Momentum and heat fluxes 



The hydrostatic pressure p = Pa jd is given by 

p = nT [l + 2 d - 2 (l + a)4> X ] , 



(28) 



where the solid volume fraction is defined by Eq. (7). To first order in the gradients, the pressure tensor Pjj is given 
by Eq. (2). As expected, while the shear viscosity n has kinetic and collisional contributions, the bulk viscosity 7 has 
only a collisional contribution. The latter coefficient is 



7 = 



22(2+1 

n{d + 2) 



^ 2 x(i + «) (i-gK 



while the shear viscosity is given by 



where 



V = Vk 



1 + + a) 



d + 2 



d + 2 r(f) 



,(d-l)/2 



(29) 



(30) 



(31) 



is the low-density value of the NS shear viscosity in the elastic limit. In Eq. (30), the subscript k denotes the 
contributions to the transport coefficients coming from the kinetic parts of the fluxes. The kinetic part rjk of the shear 
viscosity is 



nT 



Vk = 



v n ~ I Co 



l-f^(l + a)(l-3a)0X 



where the collision frequency v. q is 



^_^ £ ( 1 _ a+ | i ) (1+a) ( 1 .g). 



Here, the nominal collision frequency v is defined by 



Vo (d + 2)T{l) 
The zeroth-order contribution to the cooling rate is 



Co = ^(l- 2 )x(l + |c ), 



(32) 



(33) 



(34) 



(35) 



To first order in the spatial gradients, the heat flux q is given by Eq. (3). The thermal conductivity n and the 
coefficient u can be written as 



K — Kk 



1 + S J^2 Ml + a) 



2 2d+1 (d-l) 



(d + 2)H ^ 2 x(1 + «)(1 + ^o)ko, 
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(36) 



fj, = n k 



where 



~>d-2 



1 + 3 —2 Ml + a) 



K = 



d{d + 2) 770 
2{d- 1) m 



(37) 



(38) 



7 



is the low-density value of the thermal conductivity of an elastic gas. In addition, the kinetic parts Kk and Hk are 



d 



-j-^-K is(v K -2( ) 1 ^l + c + 3 



|i^ x (l + a) 2 [ 2 a-l + |(l + a) 



}■ 



(39) 



A*fe = — — ("k-oCo] 1 ( 1 + ^0 ln x)Co K fc + -^- c o + 3 ^ ; n , <fa(l + a) 



d(d + 2) 



x (l + Inx) [<*(<* - 1) + §(10 + 2d - 3a + 3a 2 )] J 
In Eqs. (39) and (40), Co = Co/^ 3 K £ = Kfc/«07 and the collision frequency v K is 



1 + a 



d-l 3,, „ w „ . 4 + 5d- 3(4 -d)a 
_ + _( d + 8)(l-a) + n^p^eo 



(40) 



(41) 



B. Cooling rate and fourth cumulant 

The cooling rate (, to first order in the gradients, is expressed by Eq. (4), where Co is given by Eq. (31) and 



Cv = -- d x{l-o?)[2^ 



d + 2 
128 



ci . 



The coefficient c\ is 



Cl 



4A - 2 d+1 {d + 2) X 0(1 + a) (1 - 3a) (l + f ) 
z/ c + i|d(d + 2)2 x (l- a 2 ) 



where 



(42) 



(43) 



A = 2 d - 3 <p X (l + a) {5 + 4d(l - 3a) - 9a + 3a 2 - 15a 3 - ^ [15a 3 - 3a 2 + 3(4d + 15)a - (20d + 1)] } , (44) 



= -^^X( 1 +«) [32(d- 1) + 3(1 - a)(10a 2 + 10d + 39)] . 



Finally, the fourth cumulant c is given by 

co 



32(1 - a)(l - 2a 2 



9 + 24d - a(41 - 8d) + 30(1 - a)a 2 ' 
Equation (46) coincides with the one derived by van Noije and Ernst 27 in the homogeneous cooling state. 



(45) 



(46) 



C. Dilute granular gas 



An interesting particular case corresponds to the low-density limit (</> = 0). In this case, 7 = Cu = and the NS 
transport coefficients become 

nT 



V 



(47) 



d-l 1 + c 
d v K - 2Co 



n v K 



Co 



(48) 



(49) 



Here, i/ v , Co, and v k are given by Eqs. (33), (35), and (41), respectively, with \ = 1- Equations (47)-(49) agree with 
those recently obtained 24 from the Boltzmann equation. In addition, the results reported here coincide with those 
recently derived by Kremer and Marques 26 for a dilute granular gas when one neglects the cooling effects on the 
granular temperature (i.e., when one formally takes Co = in Eqs. (47)-(49)). 
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FIG. 1: (color online) Plot of the magnitude of the first-order contribution £u to the cooling rate versus the solid volume 
fraction <j> for hard spheres (d = 3) and three different values of the coefficient of restitution a: a — 0.9, a = 0.7, and a = 0.5. 
The solid lines are the results obtained from Grad's method, while the dashed lines correspond to those obtained from the 
Chapman-Enskog method. 



V. COMPARISON WITH OTHER RESULTS 



The results derived in the preceding section provide the expressions of the NS transport coefficients and the cooling 
rate obtained by solving the inelastic Enskog equation by means of a fourteen moment method. More specifically, the 
bulk (7) and shear (77) viscosities are given by Eqs. (29) and (30), respectively, the thermal conductivity k is given by 
Eqs. (36) and (39), the coefficient \i is given by Eqs. (37) and (40), and the cooling rate £tj is given by Eqs. (42)— (45). 
It is quite apparent that the reduced forms 77*, 7*, k*, and fi* of the above transport coefficients [see Eqs. (5) and (6)] 
present a complex dependence on both the coefficient of restitution a and the solid volume fraction (f>. In order to get 
their explicit forms, the dependence of the pair correlation function x on <j) must be chosen. In the three-dimensional 
case (d = 3), a good approximation for \ is given by the Carnahan-Starling expression 29 

xi<t>) = ^ (50) 
while for hard disks [d = 2), \ can be approximated by 30 

A comparison with the results 8 obtained from the Chapman-Enskog method for inelastic hard spheres (d = 3) 
shows that the expressions of the NS transport coefficients 77, 7, k, and fi are the same as those obtained here from 
Grad's moment method. In the case of arbitrary number of dimensions d, the above expressions also agree with those 
obtained first by Lutsko 9 and more recently by the author of the present paper. 31,32 On the other hand, the first order 
contribution Cu to the cooling rate [see Eq. (42)] is different from the one derived in the Chapman-Enskog theory. 
However, a study of the dependence of (u on <f> and a shows that the results obtained from both methods are very 
similar. To illustrate these differences, Fig. 1 shows the magnitude |£u| versus the solid volume fraction (j> for three 
different values of the coefficient of restitution (a — 0.9, 0.7 and 0.5) in the case of hard spheres (d = 3). We observe 
that the Chapman-Enskog and Grad predictions for |£u| are indistinguishable in all the range of values of <f> analyzed. 
A similar behavior is also found for a two-dimensional system. Consequently we can conclude that both methods 
essentially yield the same results in the NS regime. 

Apart from comparing the present results with those derived from another different method (Chapman-Enskog 
expansion), it is quite instructive to make a comparison between the results of this paper with those previously 
reported in the literature by employing a different version of Grad's method. As mentioned in the Introduction, 
Grad's 13-moment method was already used many years ago by Jenkins and Richman to determine the NS transport 
coefficients of a dense gas of inelastic hard disks 20 and spheres. 21 Their explicit results are displayed in Appendix D 
for the sake of completeness. As usual in the conventional Grad's method, these authors did not include the fourth 
order polynomial E(V) in the trial solution (23) and so the coefficients cq and C\ vanish in their approximation. 



9 



n 




FIG. 2: (color online) Plot of the reduced shear viscosity if (a, <j)) = v( a i 0)A?(1j </>) versus the coefficient of restitution a for two 
different values of the solid volume fraction (4> = and 4> = 0.2) in the cases of hard disks (left panel) and hard spheres (right 
panel). The solid lines are the results derived here, while the dashed lines are the results obtained by Jenkins and Richman. 20 ' 21 




FIG. 3: (color online) Plot of the reduced thermal conductivity k*(cl, 4>) = n(a, 4>)/n(l, <f>) versus the coefficient of restitution 
a for two different values of the solid volume fraction (0 = and <j> = 0.2) in the cases of hard disks (left panel) and hard 
spheres (right panel). The solid lines are the results derived here, while the dashed lines are the results obtained by Jenkins 
and Richman. 20,21 




FIG. 4: (color online) Plot of the reduced coefficient fi*(a,(j)) = n/j,(a, 0)/T7v(l, (j>) versus the coefficient of restitution a for two 
different values of the solid volume fraction (0 = and <j> = 0.2) in the cases of hard disks (left panel) and hard spheres (right 
panel). The solid lines are the results derived here, while the dashed lines are the results obtained by Jenkins and Richman. 20 ' 21 
Note that the coefficient [i vanishes in the Jenkins-Richman theory for a dilute gas (0 = 0). 
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FIG. 5: (color online) Plot of the magnitude of the first-order contribution fu to the cooling rate versus the coefficient of 
restitution a for hard spheres (d = 3) and two different values of the solid volume fraction </>: (j> — 0.1 and (j> — 0.3. The solid 
lines are the results obtained here, while the dashed lines correspond to those obtained by Jenkins and Richman. 21 

A careful comparison between the results displayed in Sec. IV with those provided in Appendix D shows that the 
expressions for ij, 7, k, and [i are different in both approaches. The differences are essentially due to the assumptions 
made in Refs. 20,21 since the authors neglect the time dependence of temperature due to collisional cooling (which 
is formally equivalent to take Co = 0) and, as said before, they do not include the fourth degree polynomial E(V) 
in Grad's solution (23) (which is equivalent to take Co = Ci = 0). It is important to remark that, while the latter 
simplification is in general not relevant (except for very small values of the coefficient of restitution), the former 
assumption turns out to be quite significant beyond the quasielastic limit (i.e., for finite values of the coefficient of 
restitution). 

In order to illustrate the quantitative differences between the Jenkins- Richman theory 20,21 and the results displayed 
in Sec. IV, Figs. 2, 3 and 4 show the a-dependence of the reduced transport coefficients rf{a,4>) = r](a, (j>)/r)(l, 4>), 
K*(a,4>) = n(a, 4>)/k{1, 4>), and /j,*(a,4>) = n/j,(a, 4>)/Tk(1, <f>), respectively. Two different values of the solid volume 
fraction <j> have been considered: a dilute gas (<f> — 0) and a moderately dense gas (4> — 0.2). Moreover, for the sake 
of completeness, the above reduced transport coefficients have been plotted for disks (d = 2) and spheres (d = 3). 

It is apparent from Fig. 2 that the qualitative dependence of the shear viscosity on dissipation is relatively well 
captured by the Jenkins- Richman theory, especially in the three-dimensional case. In this case (d = 3), both Grad's 
solutions predict that while 77* is an increasing function of dissipation for dilute gases, the opposite happens at 
moderate densities. However, at a more quantitative level, there are important discrepancies between both theories, 
especially as the coefficient of restitution a decreases. The differences for the heat flux transport coefficients are much 
more significant than those observed for the shear viscosity. In the case of the thermal conductivity, Fig. 3 shows that 
both theories predict different qualitative behavior for a dilute gas (0 = 0): while the present theory shows that the 
latter coefficient increases with decreasing a, the opposite happens in the Jenkins-Richman theory. As the density 
increases, although the differences between both theories are smaller than that of a low-density gas, they are still large 
for finite dissipation, especially in the case of hard disks. As expected, the discrepancies are much more important 
for the coefficient /i (which vanishes in the elastic limit). In particular, the coefficient \x becomes negative in the 
Jenkins-Richman theory [see Eqs. (D4) and (D8)] for inelastic dense gases (</> ^ 0), but its magnitude is practically 
zero. This drawback is not present in our results since \i is always positive for any value of a and <fi. In addition, 
according to the results presented here, although the magnitude of /x is in general smaller than that of the thermal 
conductivity k, we observe that the influence of /1 on the heat transport could not be considered as negligible as the 
degree of dissipation increases. 

Finally, the dependence of the magnitude of the first order contribution to the cooling rate \(u\ on the coefficient 
of restitution is plotted in Fig. 5 for hard spheres and two values of the solid volume fraction. Note that the results 
derived by Jenkins and Richman 21 for £tj coincide with those obtained here when one neglects the coefficient c\ defined 
by Eq. (43). Figure 5 shows that the agreement between both theories is quite good, indicating that the influence 
of the coefficient c\ on |£u| is relatively small. Thus, in contrast to the behavior observed in some of the previous 
transport coefficients, the expression for (u obtained by Jenkins and Richman 21 can be considered as reliable, even 
for finite values of dissipation. 
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VI. CONCLUSIONS AND PERSPECTIVES 

In this paper, the NS transport coefficients of a granular fluid have been determined by solving the inelastic Enskog 
equation by means of Grad's moment method. As in previous works on dilute granular gases, 24-26 the solution 
proposed here differs from the conventional thirteen moment method by the inclusion of the full contracted fourth 
moment c defined by Eq. (22). To first order in the spatial gradients, this moment can be written as c = co+Ciz/ _1 V- U, 
where v is an effective collision frequency defined by Eq. (34), cq is given by Eq. (46) and the expression of c\ is given 
by Eq. (43). While the fourth cumulant Co takes into account the contributions to the transport coefficients coming 
from non-Gaussian corrections to the homogenous cooling state, the presence of c\ gives rise to a new first-order 
contribution to the cooling rate. The coefficient Co vanishes for elastic collisions while c\ = in the limits of elastic 
dense spheres (a = 1 but <j> ^ 0) and of dilute inelastic spheres (a^l but </> = 0). Thus, although both coefficients are 
different from zero for inelastic dense fluids, their impact on the NS transport coefficients is in general quite small and, 
consequently, one can neglect their contributions to the NS transport coefficients even for finite degree of dissipation. 

On the other hand, in contrast to the previous works 20,21 ' 25,26 for granular gases from Grad's moment method, 
the present results take also into account the time dependence of the granular temperature due to cooling effects. 
This dependence is accounted for by the zeroth-order cooling rate Co given by Eq. (35). In fact, if one neglects the 
small corrections due to Co and c\, the present results agree with the previous ones for dense gases 20,21 when one 
takes Co = in the expressions displayed in Sec. IV. The quantitative variation of the (reduced) transport coefficients 
on both the coefficient of restitution and density has been widely illustrated in Figs. 2-4. The comparison with the 
Jenkins- Richman theory 20,21 clearly shows that in general the influence of Co on rj* (see Fig. 2), k* (see Fig. 3), 
and /I* (see Fig. 4) is significant and so the cooling effects on granular temperature cannot be neglected beyond the 
quasielastic limit (a ~ 1). 

Furthermore, the explicit expressions obtained here for the NS transport coefficients and the cooling rate agree 
completely with those obtained several years ago from the Chapman-Enskog expansion by considering the first Sonine 
approximation. 8,9 As for ordinary gases, 5 ' 6 this agreement shows the equivalence between both approximate methods 
to solve the Enskog equation for granular dense gases in the NS regime. It must remarked that the inclusion of more 
velocity moments (for instance, all the third and fourth degree velocity moments) would change the final results since 
for instance there would be likely additional contributions to the cooling rate. Nevertheless, for practical purposes, 
the inclusion of those new terms in Grad's solution (23) makes analytic calculations much ore difficult since higher 
order collision integrals should be evaluated to compute the new contributions to the momentum and heat fluxes. 

The derivation of explicit expressions for the transport coefficients is perhaps one of the most important challenges 
of granular gas research. The theoretical results reported in this paper cover part of this challenge, at least in the NS 
domain. Nevertheless, the present results have some restrictions. First, although the Enskog equation retains spatial 
correlations arising from volume exclusion effects, it still assumes uncorrelated particle velocities (molecular chaos 
hypothesis). Therefore, it is expected that the results reported here only apply to moderate densities (solid volume 
fraction typically smaller than or equal to 0.25). However, in spite of this limitation, there is substantial evidence 
in the literature 33 on the reliability of the Enskog kinetic theory to accurately describe macroscopic properties (such 
as transport coefficients) for a wide range of densities and/or collisional dissipation. Another important limitation is 
the accuracy of the results obtained here for quite extreme values of dissipation. As mentioned in the Introduction, 
although the approximate expressions for the NS transport coefficients displayed in Sec. IV compare in general quite 
well with computer simulations, 13 there are significant discrepancies for the heat flux transport coefficients for small 
values of a (say, for instance, a < 0.7). Therefore, the reliability of the expressions (36) and (37) for the thermal 
conductivity k and the coefficient fi, respectively, can be questionable for this range of small values of the coefficient of 
restitution. In the context of the Chapman-Enskog solution, the above discrepancies between theory and simulation 
can be in part mitigated 16,34 when one uses the homogeneous cooling state distribution instead of the Maxwcllian 
distribution as the weight function in the corresponding first-Sonine approximation. The modification of Grad's 
method employed here by replacing the Gaussian distribution /m in the trial solution (23) by the homogeneous 
cooling state distribution is an interesting open problem to be studied in the near future. This will allow us to offer 
a Grad's solution that will be reliable even for extreme values of dissipation. 
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Appendix A: Collisional contributions to the fluxes 



The collisional transfer contributions to the pressure tensor and heat flux are determined from Eqs. (19) and (20), 
respectively. Consider first the pressure tensor which becomes up to first order in the gradients 



pc 1 + a d 

Pij = ^ TfL(7 



xjdVij dV 2 J dSQ(d- -g)(5 ■g) 2 a i a J [/o(V 1 )/o(V 2 ) + / 1 (V 1 )/ (V 2 ) 



+/o(Vi)/i(V 2 ) - ^/o(V 2 )5 • V/oCVO + ^/ (Vi)5 • V/ (V 2 ) 



where 



/o(V)=/ M (V) (l + ^(F)), 

is the part of the trial distribution (23) of zeroth-order in spatial gradients and 

2 to fc Ci 



/i(V) = / M (V) 



— viVjiUj + d + 2nT2 



(Al) 



(A2) 



(A3) 



is the part of the trial distribution (23) of first-order in spatial gradients. Upon writing Eqs. (A1)-(A3) use has been 
made of Eq. (27). The integration over solid angle in Eq. (Al) yields 



nc 1 + a B 2 a 

P ii = — a i — — to<t 

13 4 d + 2 



dUtl + a B 3 d+1 
where the coefficients B n arc defined by 



X / dVi I dV 2 [/o(V 1 )/o(V 2 ) + 2/ 1 (Vi)/ (V 2 )] (2g i9j + g 2 S tJ ) 

J dVi y dVa/oCVO ^^ g- 1 [&g jgfc + g 2 (g fc % + gi g jfc + gj S ik )\ , (A4) 



B„ = y da6(?-g)(a-g)"=7r( d - 1 



)/2 i I 2 / 



(A5) 



The expression (A4) can be more explicitly written when one takes into account the forms of /o and /i. After some 
algebra, the result is 



1 + a 



na X nT 



B 9 



d+2 



(1 + a)na d X ^ 



0U e l + a B 3 d+1 



J dVi J rfV 2 /o(Vi)/ (V 2 )^- {g 1 [gigjgk + g 2 {gAj + 9 iS 3k + gjS lk )} } 



= B 2 ^^na d X nT + + a)na d X n i3 



B, 



2 

d+1 



-ma 



d + 2 

d+1 1 dVi I wtfotyMoiv^ 



4d 2 
1 + a 
2~" /v '" ' d + 2 
d+lT(^ 



^Uui + diUj-^V-V) -+X-V: 



B 2 ^^na d X nT + -^(1 + a)na d X % 



^Uui + diUj-lSijV-TJ) +S ii Y-V 



• (A6) 



Note that in the last term of Eq. (A6) nonlinear terms in Co have been neglected. From Eq. (A6) one can easily identify 
the collisional transfer contributions to the hydrostatic pressure p c , the shear viscosity r] c and the bulk viscosity 7. 
They are given by 



Pc = 2 d - 2 {l + a)cj )X nT, 



(A7) 



02(<J-l)r (d\ , 



(AS) 
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2 d_1 d 
Vc = TT-^X(1 + + -r—rj, 



d+2 



d+2 



(A9) 



where use has been made of the definition (7) of the solid volume fraction (j). 

The collisional transfer contribution to the heat flux to first order in the gradients can be obtained in a similar way. 
The result is 



<n = 



l + a 
2 

3 B 2 



ma d X J dVi j dV 2 J da 0(5 ■ g)(5 ■ g) 2 {a- G)a, 
B 3 m<7 d+1 



/i(Vi)/ (V 2 ) + ^/o(Vi)ff • V/ (V 2 ) 



2d + 2 

3_B2_ 
2d+2 



,Wy./; - •■ /- —X(1+«) j dV x J dV 2 /o(Vi)/ (V s 



s 3 r(^) /_ . 7 \ o j 1 1 IT 



2d r(|) V 1 + ^ C0 J n2(Td+1 V^ (1 + a) - 



.9- 1 (g.G) 2 + 5 G 2 + |g(g.G) + i g 3 



(A10) 



From Eq. (A10) one may identify the collisional contributions to the thermal conductivity k c and the coefficient [i c . 
They are given by 



2 d ~ 2 2 2d_3 r (-) ( 7 



(All) 



d+2 



(j)X(l + a)fi k . 



(A12) 



Appendix B: Kinetic contributions to the fluxes 

The kinetic contributions to the momentum and heat fluxes are defined by Eq. (18). To obtain them in the NS 
approximation, one has first to expand the collision operator Je[/, /] to first order in the gradients. To do so the 
following results are needed: 



1 3 In x 

X(r,r ± cr\n) ->• % ( 1 ± 2 n ~g^~ (T ' Vlnn 



(Bl) 



/(r ± a, V, t) -> / (r, V, t) + h (r, V, t) ± a ■ V/„(r, V, t), 



(B2) 



where /o and /i are defined in Eqs. (A2) and (A3), respectively, and \ is obtained from the functional x( r , r ± 
by evaluating all density fields at n(r, t). The collision operator to first order then becomes 8 



■h'j.l] - i • io^Ajx;./:: -vinn • i/c 



6> 



5/o 



5/ 



dV,, 



8V 
V-U-£/i, 



(V/o) 



■ VlnT 



where 



and we have introduced the operators 



(0) 



and 



j<°> [a, y] = x^- 1 y dv 2 y rf?e(? • g )(5 • g) [a- 2 AK)y(v 2 ) - x( Vl )Y(v 2 

/C[A] - A y rfv 2 y ^6(? • g)(? • g)<r [ a - 2 /(°)(v' 1 )A(v 2 ) + /( )( Vl )A(v 2 ) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 
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Here, and v 2 are denned by Eq. (11). 

We consider the kinetic contributions to the shear viscosity 77. Multiply both sides of Eq. (8) by mViVj and integrate 
over velocity to get 

dtP* + f*V ■ U + UtdtF* + PfjdtUi + P&eUj + jL^dtitfSji + q$6u + q\k 3 ) = J dvmV^Mf, /], (B8) 



where use has been made of the result 



J dv mViVjVefiV) = + QjSu + 



The integral on the right hand side of Eq. (B8) can performed by using the definition (B3) 

'dfo 



J dv mViVjMf, /] = - J dv mViVjCh + \J dv rn\^V 3 IC e 



dV n 



dv mViVjJCi 



dfo 



(B9) 



V-U. (BIO) 



The collision integrals appearing in Eq. (BIO) can be evaluated by considering the explicit forms of f Q and f\. The 
first integral is given by 8 ' 24 ' 31 



J dv mViVjCfx = Vr,Uij + nT<Co%, 



(Bll) 



where v v and Co are given by Eqs. (33) and (35), respectively. Note that nonlinear terms in Cq, Ily and q fc have been 
neglected in Eq. (Bll) when f and f\ are replaced by their Grad's approximations (A2) and (A3), respectively. The 
remaining two integrals can be performed using the definition of K in Eq. (B7). The first integral is 



/ 



dv mViVjlCi 



dfo 



dV v 



A ip = g\ J dVi J dv 2 mv li v lj J d*e(* ■ S ){a ■ g )d t 



dV 2 



2p 



A fr> . 



(B12) 



A simpler form of this integral is obtained by changing variables to integrate over and V 2 instead of Vi and V 2 
in the first term of Eq. (B12). The Jacobian of the transformation is a and <x • g = — aa ■ g'. Also, Vi(V' 1 , V 2 ) = 
V" = Vi — |(1 + a)er(er ■ g). The integral then becomes 



dvmViVjJCi 



dfo 



dV v 



Hp 



m X a d J dV 1 J dV 2 .h(V 1 ) df ^ 2) J dZQ{Z ■ g)(S • g)a £ (V{X' 3 ~ V u V l3 ) A, 



dV 2p 

= m X ^* d J dVi J dV 2 fo(V 1 ) df{ ° d ) ^ 2) J d*Q(* ■ • g 
x [2 (Vudjdt + Vijdidi) - (1 + a){a ■ g)<7i<7,-3*] A £p 

' f V J dVi j dV 2 /o(Vi)/ (V 2 ) J d*e(* ■ g)(5 • g) 



x [4 {Vndjdfdp + Vijdidedp) - 3(1 + a) (a ■ g)didjdidp] Ae p 
2 d-i 



d + 2 



nT0x(l + a)(l-3a)A«, 



(B13) 



where in the last step use has been made of the angular integrals 



/ 



da@{a • g) (ct • g,)aidjdt 



B 2 
d + 2 



(Sijge + 5u9j + S 3 e9z) , 



(B14) 



B 2 



[2 (gigjSep + gtgi5 pj + g l g p 5i j + gegjS ip + g p g ] 5 l i> + gig p 5 l 



{d + 2)(d + 4) 

+g 2 (Sijdfp + 5uSjp + S ip 5ji)] . (B15) 
The last integral appearing in Eq. (BIO) can be performed by using similar mathematical steps. It is given by 

-0/(0)" 



j dv mViVjJCi 



dVe 



6 l3 2 d - 2 nTx(f>(l + a)(l - 3a). 



(B16) 
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We are interested in the solution to Eq. (B8) in the NS approximation. Thus, in order to solve (B8), one needs to 
make use of the balance equations (15)— (17) up to first order in the spatial gradients: 

d t n^ -U-Vn-nV-U, d t U -> -U • VU - p _1 Vf>, d t T -> -U • VT — r~pV • U — (T, (B17) 

an 

where the pressure p is given by Eq. (28) and the cooling rate can be written as (see Appendix C) 

C^Co + (Cio + CiiCi)V-U. (B18) 

According to Eq. (2), the kinetic contribution rjk to the shear viscosity is defined as 

Iltf = -fJfcAy. (B19) 

The coefficient r/k can be easily obtained from Eq. (B8) when one takes into account Eqs. (BIO), (Bll), (B13), and 
(B16). The corresponding equation for r/k is 



(dt + v v )Vk =nT 



1_ Jr^ (1+a)(1 ~ 3a) 



(B20) 



where the time derivative dtrjk must be evaluated at zeroth-order in spatial gradients. From dimensional analysis 
r/k oc T 1 / 2 and so, d t r\k — \r\k&t InT = — ^(oVk- With this result, the solution to Eq. (B20) is given by Eq. (32). 

Apart from obtaining the coefficients in Eq. (B8) proportional to the divergence of the flow velocity allows one 
to determine Cio- It is given by 



Cio = -3- 



d 



(B21) 



This first-order contribution to the cooling rate will be also determined in the Appendix C by following a different 
route. Moreover, the coefficient c\ = at this level of approximation. A nonzero contribution to c\ will be obtained 
when we determine the (contracted) fourth degree velocity moment of /. 

The kinetic parts of the thermal conductivity k and the coefficient fi are obtained in a similar way. Multiplication 
of Eq. (8) by '^V 2 Vi and integration over the velocity yields 



dtq" 



2 

d + 2 



nTd t U l + 7^(nT 2 ) 
Zm 



d+2 
2m 



dj (nT 2 ) 



d + 2 



m , 



nTU-V[/,= / dv-V 2 ViJ E [f,f] 



(B22) 



where only linear terms in the spatial gradients have been considered on the left hand side of Eq. (B22). In addition, 
I have assumed that c is uniform (see Eq. (46)) and have used the relation 



/ 



d^V 2 V l .f{Y)= i 



d + 2nT^ 



Sij (c + Cl is X V-U) + — 



T fd + A 



(B23) 



m m 

Equation (B22) can be more explicitly written when one takes into account the balance equations (B17). The result 
is 

d + 2T 2 



m 



1 , Co * a 9 P* 



) d t n + (2 + co - P*) d t T = J dv^V 2 V t ML f], (B24) 



where 



nT 



= l + 2 d - 2 4> X (l + a). 



(B25) 



The collision integral on the right side of (B24) can be evaluated by using similar mathematical steps as those made 
before in Eqs. (B11)-(B16). After a tedious and long algebra, one gets 8,24 ' 31 
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dj In n 
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+ -d j \nT 



d v'J± V 2 Vi Kj 
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nT 2 
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dilnT, (B26) 
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where v K is given by Eq. (41). As before, linear terms in Co, Ily, q fc and c\ have been only retained in Eq. (B26). 
The constitutive equation for q fc is 



(B27) 



Dimensional analysis shows that Kk oc T 1 / 2 and u/j cx T 3 / 2 . Thus, according to Eq. (B27), to first order in the 
gradients the time derivative of the kinetic contribution to the heat flux can be written as 



dtqf = 2Co«fcftT + Co 



Tkh f dhi X \, 3 



where use has been made of the relation 



di(d t T) -+ -di(CoT) = -CoT ( 1 + 4>^- J ft Inn - -CoTft InT. 



(B28) 



(B29) 



Substitution of Eqs. (B26) and (B28) into Eq. (B24) leads to the following set of equations when one equates coefficients 
pertaining to the density and temperature gradients: 



d -\-2 riT nT 

K^k + -z (2 + c - V*) = v^k - 2 d ~Vx(i + a) — 

2 m m 



2{d + 2) + 3(1 + a)(2a - 1) + -c (l + af 
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d In y 
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d + 2T 2 f c , 



2(d + 2) + 3a(a - 1) + ^(10 + 2d - 3a + 3a 2 ) 



(B31) 



The solution to Eqs. (B30) and (B31) gives the expressions (39) and (40) for Kk and /ik, respectively, once the 
explicit form (B25) for p* is considered. 



Appendix C: Cooling rate £ and fourth cumulant Co 



The cooling rate C is defined by Eq. (21). Up to the first order in gradients, ( is 
1-a 2 



where 



c = 



AdnT 



ma d -' X \ dVi / rfV 2 / d5e(S-g)(S-g) d /o(Vi)/ (V 2 ) 
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Cio^-^^^m^-^l dV! | dV 2 ff 2 /o(V!)/ (V 2 ), 



(CI) 
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(C3) 



Cn = -B3 



1-Q 2 



m^-V" 1 / dVi / dV 2 5 3 S(V 2 )/o(Vi)/o(V 2 ). 



(C4) 



Here, the coefficients -B„ are defined by Eq. (A5) and v is given by Eq. (34). The integration over velocity can be 
carried out in Eqs. (C2)-(C4) when one takes into account the expression (A2) of /o- Neglecting nonlinear terms in 
cq, the coefficients Co and C10 are given by Eqs. (35) and (B21), respectively while Cn is 



3(rf + 2) 2 A ,3cd\ 

Cll = ^r (1 - a)x 1+ 64 '• 



(C5) 
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The complete determination of the cooling rate and the NS transport coefficients still requires to get the quantities 
Ci and Co- To evaluate them, one multiplies both sides of Eq. (8) by V 4 and integrates over velocity. In the NS order, 
one gets 



d(d + 2) 



(i + f + £v.u) 



d t {nT 2 ) + ^±inT 2 V • U + U • V(nT 2 ) 



J dvV 4 J E [f,f}. (C6) 



The left hand side of Eq. (C6) can be simplified when one takes into account the balance equations (B17). To first 
order in spatial gradients, the result is 



nT 2 

-2<d + 2) — 



Co+(^(p*-l) + Cio + CiiCi)V-U 



^ Cl v-uJ = J dvV 4 M.f,.f}. 



(C7) 



The collision integral on the right side can be computed by using similar mathematical steps as those made before. 
Neglecting nonlinear terms in cq and C\ and after a tedious algebra, one gets 



/ 



d+-+a A 



+ 



rrr 



V-U, 



(C8) 



where A and v$ are given by Eqs. (44) and (45), respectively. Upon deriving Eq. (C8), use has been made of the 
partial results 



/ 



dw V A Ch = v ( 



nT 2 



(C9) 
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dv V A lCi 



dX 



nT 2 
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Substitution of Eq. (C8) into Eq. (C7) allows one to explicitly determine Co and c\. They are given by Eqs. (43) 
and (46), respectively. Note that in Eqs. (43) and (46) nonlinear terms in c and C\ (like c 2 , c a ci and c\) have been 
neglected. 

Appendix D: Jenkins and Richman expressions 

The expressions obtained years ago by Jenkins and Richman 20,21 for the NS transport coefficients of a dense gas of 
inelastic hard diks 20 and spheres 21 are displayed in this Appendix. These authors solved the Enskog kinetic equation 
from thirteen Grad's moment method. 

In the case of smooth hard disks (d — 2), their results are 20 
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Here, 7/o and kq are given by Eqs. (31) and (38), respectively. 
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In the case of smooth hard spheres (d = 3), their results are 21 
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